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A detailed study of the high-mass clump interacting with 
the bubble N10 
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Abstract We performed a detailed study of the high- 
mass clump interacting with bubble N10 based on the 
spectral lines 12 CO(3-2), HCO + (4-3), N 2 H + (4-3) 
and CH3OH(7(0, 7) — 6(0,6)) and continuum emission 
data at 450 /^m and 850 /im released on CADC and 
Spitzer data. Blue-shifted optically thick line 12 CO(3 — 
2) seems to indicate that the outer envelope of the high- 
mass clump is still falling toward the center. Detection 
of C*iJ 3 OiJ(7(0, 7) - 6(0,6)) suggests that a hot core 
has formed around YSO N10-7. And position-velocity 
diagram of iV2-ff + (4 — 3) indicates the cold dense core of 
the clump has not been destroyed by the star formation 
activities. The mass of N10-7 is about 27.44 M Q . The 
ratio HCO + (4 - 3)/A r 2 iJ+(4 - 3) in the outer part of 
the clump is larger than that in the inner part of it. 
The reason may be that the CO abundance relative to 
A^2-ff + (4 — 3) increased in the outer part of the high- 
mass clump, more A^2-ff + (4 — 3) were converted into 
HCO+(4-3). 
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1 Introduction 

The expansion of the HII regions is extreme interest to 
studies of star formation as their expansion may trigger 
new generations of star formation into being within the 
molecular mater ial surrounding the bubbles (Thomp- 
son et al. | |2012l ) ). Evidence of triggering has b een re - 
ported by many autho rs (e.g. Deharv eng et al. | (|2010h : 
Zavagno et al. | |2010l ); Kang et al. | |2009l) ). It should 
be noted, the majority of observational studies into trig- 
gered star formation near SNR or HII regions take a 
phcnomenological approach, the evidence of triggered 
star fo rmation is not very conclusive (Kendrew et al. | 
( 20121) ). The statistical approach may address the un- 
certainties inherent in observations of individual HII re- 
gions. One detailed statistical study of massive star for- 
mation in the environment of 322 Spitzer mid-infrared 
bubbles by using the Red MSX source survey for mas- 
sive young stellar objects (YSOs) suggest that the frac- 
tion of massive stars in the Milky Way formed by trig- 
gering could be be tween 14 and 30 per cent ( Thomp- 
son et al. | |2012l) ). Kendrew et al. J |2012l )) made 
a similar statistical study with 5106 infrared bubbles, 
they estimated that approximately 22 per cent of mas- 
sive young stellar stars may have formed as a result of 
feedback from expanding HII regions. Therefore, the 
infrared dust bubbles could be good sites for us to find 
high-mass YSOs and study the process of high-mass 
star formation. 

N10 is a bright mid-infrared and radio continuum 
bubble with an elliptical or slight cometary shape 
with a n opening at southeast of bubble (Watson et alj 
( 20081) ). It s kinem atic distance is about 4.6 kpc (Dehar- 
veng et all (|2010h ). The bubble is bordered on two sides 
by infrared dark clouds (IRDCs), they are in teracting 
with the HII regions (Deharveng et al| |2010h ). Espe- 
cially for the IRDC in which one medi um-to -high mass 
YSO has been found by Watson et al| |2008l) . the YSO 
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is coinc ident with one very dense dust core (Miettinen 
| (|2012l )). So this IRDC is one high- mass clump and 
is forming YSOs. It provide us a opportunity to study 
initial conditions and processes of YSO's formation and 
triggered star formation. 

In this paper, wc study the high-mass clump and 
high-mass YSO formed in it based on archive data, 
which including continuum data at wavclengthcs 1100 
/im, 850 /xm, 450 /im and molecular lines of 12 CO(3— 2), 
J ffCO+(4-3), N 2 H + {4-3), aff 3 O-ff(7(0,7)-6(0,6)). 
The observational details and data analysis are pre- 
sented in the following sections. 

2 Data 

Three large-scale survey results were used in our 
study, which includes Galactic Legacy Infrared Mid- 
Plane Survey Extraordinaire (GLIMPSE), the NRAO 
VLA Sky Survey (NVSS) and BOLOCAM Galactic 
Plane Survey (BGPS). GLIMPSE is performed using 
the Spitzer Space Telescope by Spitzcr-Infrarcd Array 
Camera (IRAC). It is a mid-infrared survey of the inner 
Galaxy at 3.6, 4.5, 5.8 and 8 /im with angular resolu- 
tion of between 1'.'5 and l'/9, and a very large catalog 
of star s wer e extracted from the survey (Churchwcll 
et al| |2009l) ). Here we used the mosaicked images 
from GIMPSE acquired by Spitzer-IRAC at 8 /im and 
the Point-Source catalog. NVSS is a 1.4 GHz radio 
continuum survey covering the entire sky north of - 
40 deg declination with a n angular resolution of about 
45'.'(Condon et al| |l998h ). BGPS is a 1.1 mm contin- 
uum survey of the Galactic Plane made using Bolocam 
on the Caltech submillimeter Observatory, the survey 
coverage totals 170 square degrees (with 33'.'FWHM 
effective resolution) and detecte d about 8400 point 
sources (Rosolowsky et al| |2010[) ). 

In addition, we used the James Clerk Maxwell Tele- 
scope (JCMT) data released on Canadian Astronomy 
Data Centre (CADC), they including 12 CO(3 - 2), 
HCO+(4: - 3), N 2 H+(4: - 3), CH 3 OH(7{0, 7) - 6(0, 6)) 
data and dust continuum emission at 450 fxui and 850 
/im. The co-added spectral cubes were binned to 0.2 
fcms -1 , and the data were smoothed with 6'.'Gaussian 
kernel. The fin al spa tial resolution of each cube is 
16'.'(Smith et al| d2008h ). The angular resolution of 450 
/itm and 850 jim is 8'.'5 and 14'.'5 respectively (Leeuw& 
Robsonl |2009h L 
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Fig. 1 N10, the 12 CO (3-2) contour levels range from 
90.0 to 159.5 by 9.93 Kkms -1 are superimposed on the 
composite colour image with 20 cm emission in red and 8 
/im emission in green. The high-mass clump which we study 
in detail was marked with one red square. The whi te sta r 
denotes the YSO (N10-7) identified by Watson et al. | (120081 ). 

3 Results and discussion 

3.1 N10 

Figure 1 shows the 12 CO(3 — 2) emission contours over 
the two colour image with 20 cm emission in red and 
8 fj,m emission in green. The border of N10 is well 
traced by the 8 ^m emission. The 20cm emission in 
the bubble indicates that there is a large HII region 
ionized by central stars. The high-mass clump, which 
we will study in detail, was marked with one red square. 
An obvious arch border appears between the bubble 
and the high-mass clump, which suggests that there 
exists strong interaction. One very compact CO cloud 
is coincident with the high-mass clump, the medium- 
to- hig h-mass YSO (N10-7) identified by Watson et al. | 
(|2008l ) is located at the center of it. The more extended 
CO cloud appear on the upper left to the center of N10. 
The CO emission toward the center of the bubble in 
projection could be arising from the front or back side 
of the bubble. 

We estimated the dynamical age of N10 us ing the 
model described by Dyson & Williams | |l980t ) with a 
given radius R as: 

t(R) = ^ - ll (i) 

7c s [ R s 

where c s is the sound velocity in the ionized gas 
(c s = 10km s" 1 ) and R s is the radius of the Stromgrcn 
sphere given by R s — (3Yi yc /47rn,Q ap)z , where Nl V c 
is the number of ionizing photons emitted by the star 
per second, n is the original ambient density, and 
ap = 2.6 x 10~ 13 cm 3 s _1 is the hydrogen recombination 
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coefficient to all levels above the ground level. Here 
we adopt a Lyman continuum photon flux of about 
2.2 x lO^phs- 1 (Watson et alj (|2008l ,) and a radius 
of 1.61 pc. For that no CO cores look like swept by 
expanding HII region on the border of N10, it is dif- 
ficult for us to estimate the original ambient density 
of N10. Massive star formations usually take place in 
the giant molecular clouds and the densities of the giant 
molec ular clouds are about 10 3 to 10 4 cm~ 3 (Goldsmithj 
1 19871 ) ). So it seems to be safe for us to assume the 
low value of original ambient density of N10 is about 
10 3 cm -3 and obtain a lower bound on the dynamical 
age, about 9.17 x 10 4 yr. Taking into consider that stars 
are formed in the dense clumps of the giant molecular 
cloud, the true original ambient density may be larger 
than that we assumed, the true dynamical age of N10 
should be larger than 9.17 x 10 4 yr. 
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3.2 The High-mass Clump 
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Fig. 2 The white solid contours denote 850 /im emission, 
contour levels range from 0.77 to 3.43 by 0.53 Jy / beam. 
The white dotted contours denote 450 (im, contour levels 
range from 3.57 to 11.95 by 2.10 Jy / beam. The lowest 
levels of 850 /jm and 450 /im are both larger than 3a. The 
background is 8 /im emission. The star represents N10-7. 

Past continuum observations at diffcren wavelength 
have proven that there is one dust condensation in the 
high-mass cl ump, the YSO (N10-7) is forming in it 
(Miettinen| (j2012f )). The JCMT continuum observa- 
tion results at 450 fim (white dotted contours) and 850 
/zm (white solid contours) display a elongated dust core 
from north to south (see Fig. 2). N10-7 is nearly at the 
center of it. No mid-IR diffuse emission and lack of 
20 cm continuum emission suggest that high-mass star 
formation is in its very early stage, or the high-mass 
YSO was embedded in the clump very deeply. 

The total dust mass of the clump is calculated using 
M = S V D 2 /k v B v (Td), where S v is the flux at frequency 




Fig. 3 12 CO(3 - 2), HCO+(4: - 3), N 2 H+(<i - 3) and 
CH-sOH(7(0,7) — 6(0,6)) contours map of the clump are 
superimposed on the 8 /im emission. The contour levels 
range from 90.0 to 159.5 by 9.93 K /ferns" 1 for 12 CO(3-2), 
3.0 to 15.78 by 1.83 Kkms- 1 for HCO + , 3.0 to 14.43 by 
1.63 Kkms' 1 for N 2 H+ and 3.0 to 17.15 by 2.02 K km s _1 
for CH 3 OH(7(0, 7) - 6(0, 6)). The star represents N10-7. 



v, D is the distance (4.6 kpc), B v (Td) is the plank func- 
tion, and k v is the dust opacity per unit gas /dust mass. 
Td is the dust temperature, it was obtained by SED fit- 
ting with a simple greybody model. The total flux of 
the clump at 450 /im, 850 fxxa and 1.1 mm are about 
99.33, 13.30 and 7.6 Jy, respectively. The flux at 450 
and 850 /im were derived from JCMT data by oursclf, 
the flux at 1.1 mm were obtained from BGPS point 
source catalog and multiplied by a factor of 1.5. The 
best fitted result indicates that Td = 27 K. Using the 
dust opacity of /cssn = 0.0 2cm 2 q~ 1 obtained from Os- 



senkopf & Henning (|1994l ). the ratio of dust to gas is 
taken as 0.01. We derived a total mass of M = 908M Q 
from 850/im flux. The area of the clump is 1333 arcsec 
square. Assuming the clump is spherical, we get the 
core radius of 20'.'6 (corresponding to 0.46 pc). Then 
the average numb er den sity is 4.59 x 10 4 cm~ 3 . Fol low- 
mg Beuther et al| (|2012h and Longmore et al| |2012l ). we 
2 and calculate Virial mass using the fo r- 
Richardson and Wolfendal«| (jl98ct . 
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where AV is the full width at half-maximum intensity 
in /cms -1 , R is the radius of the molecular cloud in par- 
sees, M is cloud mass in unit of solar mass. Since star 
is formed in the cold dense core traced by -A^-ff -1- , we 
used the Gaussian line width of N2H + obtained at the 
position of YSO (N10-7) to calculate the Virial mass 
of the clump. Finally, we obtained M V i r =1640M Q . 
This indicate the clu mp m ay be gravitational stable. 
However, Miettinen | (|2012 ) obtained a total mass of 
3335+704M Q for the clump based on 870 /im obser- 
vation with Apex submillimeter telescope. The mass 
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difference may be due to the observation sensitivity 
difference b etween the two telescopes. If the result of 
Miettinen is true, the clump tends be under 

gravitational unstable. Taking into that the clump is 
forming one YSO in it, it seems to be more reason- 
able. Just as that pointed out by Beuther et al| (|2012 ). 
there is large uncertainties for the mass estimated from 
the dust continuum emission and Virial mass calculated 
with A^2-ff + , such a result are not very conclusive. 

3.3 The kinematics of the clump 




Fig. 4 The molecular profiles extracted from the positions 
around the YSO N10-7 are plotted, the panel marked by 
a star correspond to the YSO N10-7. From up to down 
in each panel is 12 CO(3 - 2), #CO+(4 - 3), Y 2 #+(4 - 
3), CH 3 OH(7(0,7) - 6(0,6)), respectively. The velocity is 
between 25 to 50 km/s. CH 3 OH, N 2 H + , HCO + profiles 
are scaled up by a factor 5, 3 and 5, respectively, and their 
base line is added 4, 8, 12 Kkms -1 , respectively. The 
vertical line indicates the system velocity of the high-mass 
clump which is 35.64 ferns -1 . 

N2H + is not strongly affected by freeze out on grain 
surfaces, it is a good tracer of younger cold and mas- 
sive cloud cores. HCO + is a highly abundant molecule, 
with abundance especially enhanced around regions of 
higher fractional ionization, it is also enhanced by the 
presence of outflows where shock-ge nerat ed radiation 
fields are present (Vasyunina et alj ( 201ll )). CH 3 OH 
molecules arc formed in the icy mantles of interstellar 
dust grains by hydrogenation of CO molecules at tem- 
perature of 10K, as the young protostellar object evolve 
and warm up its environment the CH 3 OH molecule 
sublimate from the dust grains into the gas phase at 
100K (TorstenssorJ |201lh ). Detection of HCO+{A - 3) 
and CH3OH toward the high-mass clump indicates 
that one hot core has formed around the YSO N10- 
7. On the other hand, detection of N2H + (A — 3) in- 
dicates star formation in the clump is still in its early 



stage, the natal molecular cloud was not destroyed by 
the feedback of YSO (see Fig.3). 

Fig.4 displays the 12 CO(3 - 2), HCO+{A - 3), 
N 2 H+{A - 3), CH 3 OH(7(Q,7) - 6(0,6)) profiles ac- 
cording to their position around the YSO N10-7. Here 
CH 3 OH(7{0, 7) -6(0, 6)), Y 2 i?+(4-3) and HCO+{A- 
3) profiles are scaled up by a factor of 5, 3 and 5, re- 
spectively. The panel marked by a star indicates the 
position of YSO N10-7. the vertical line in the map in- 
dicates the system velocity (35.64 ferns -1 ) of the high- 
mass clump. We noted that the profiles of optically 
thick line 12 CO(3 — 2) include two components at veloc- 
ity 34.57 and 41.77 km s _1 . The first is the stronger one 
and is usually blue-shifted with respect to the system 
velocity. This component traces the high-mass clump, 
the blue-shifted profiles of it may indicate that the outer 
envelope traced by CO is still falling toward the center. 
The profiles of HCO + (A — 3) in the all upper two rows 
and the left two panels of the third and fourth rows are 
red-shifted with respect to the system velocity. The 
Y2-ff + (4 — 3) profile at the position of YSO show evi- 
dence of self-absorption dip. The profiles in upper two 
rows of panels are red-shifted with respect to the system 
velocity, while in the lower three rows of panels they 
become blue-shifted. The position-velocity diagram of 
Y2-ff + (4 — 3) shows a velocity gradient along the major 
axis of the high-mass clump (see Fig. 5). Such a velocity 
gradient may be a combination of both rotation and i n- 
fall motions of the cold dense core (Tobin et al| (|2012 )). 
This supports our idea that cold dense core of the clump 
has not been destroyed by the feedback of the YSO 
N10-7 yet. The optically thin CH 3 OH(7(0, 7) - 6(0, 6)) 
profiles have only one peak with obvious line wings, but 
corresponding position- velocity diagram (see Fig. 6) dis- 
play only slightly red-shifted and blue-shifted lobes. In 
some extent, this supports the exist of outflow. How- 
ever, other motions such as rotation, infall and turbu- 
lence could also cause this effect, we can not exclude 
these motions. More high-resolution observations and 
studies are needed to draw a conclusive result. 

Assuming that line wings of CH 3 OH(7(0, 7) — 
6(0, 6)) indicate outflows truly, we plot the blue and red 
lobe of the outflow (see Fig. 7), the distribution of red 
lobe is much larger than that of blue lobe. The back- 
ground is the ratio HCO+(A - 3)/N 2 H+(A - 3). It is 
interesting to find that ratio HCO+ (4-3) /N 2 H+(A-3) 
in the outer part is larger than that in the inner part of 
the clump. For that HCO + is created by the reaction 
between N2H + and CO. The reason may be that the 
CO abundance relative to N2H + (A — 3) increased in the 
outer part of the high-mass clump, more N 2 H + (A — 3) 
were converted into HCO + (A — 3). Increasing of the 
relative abundance of HCO + (4 — 3) increases the ratio 
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VELOCITY (Km/s) 



Fig. 5 PV diagram of A r 2-ff + (4 — 3) emission along the 
major axis of the high-mass clump. The contours are from 
10% to 90% with steps of 10% of the peak of the integrated 
intensity. 



of HCO+{4 - 3)/N 2 H+{4 - 3). The distribution of 
HCO + (A — 3) is nearly the same as that of red lobe. 
This is consistent with the conclusion that HCO + (4—3) 
abundance will be enhanced by the presence of outflows. 
No 20cm continuum emission was detected to the clump 
suggests that it could not be highly fractional ionized. 
Enhanced HCO+(4-3) abundance is likely caused by 
the outflow. This is also consistent with the fact that 
most profiles of HCO + (A — 3) are red-shifted with re- 
spect to the system velocity. 
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Fig. 6 PV diagram for the CH S OH(7{0, 7) -6(0, 6)) emis- 
sion along the major axis of the high-mass clump. The con- 
tours are from 10% to 90% with steps of 10% of the peak 
of the integrated intensity. The vertical black lines indicate 
the blue- and red-shifted components' velocity boundary. 



Fig. 7 The components of the CH3OH outflow su- 
perimposed on the integrated intensity ratio HCO + (4: — 
3)/JV2-H" + (4— 3) (red background). The velocity components 
is integrated in (27,32.5) and (37.4,42.5) km s" 1 for its blue 
and redshifted lobes, the blue and red contours represent 
blue lobe and red lobe, respectively. The white contours 
denote HCO + (4 — 3) molecule's emission. 



3.4 SED fitting of the central massive star 



Watson et al| (|2008h fitted the SED of the young stellar 
object (N10-7) formed in the high- mass cl ump using the 
tool developed by Robitaille et al| |2006l ) based on the 
Spitzer IRAC data, they found that N10-7 is one Stage 
I YSO with a mass of 12.4 M . However, their SED fit- 
ting results may not be very accurate due to the lack of 
data at longer wavelengthes. We fit the SED of N10-7 
using more data, which including Spitzer-IRAC data at 
3.6, 4.5, 5.8 and 8 fxm, SCUBA data at 450 and 850 fxm 
and BOLOCAM data at 1.1 mm. The visual extinction 
is estimated from the classical relations Nh+h 2 /E(B — 
V) = 5.8 x VS 21[ part%clescTn~ 2 mag~ 1 (Bohlin et al| 



(|!978h ) and A v = 3.1E(B - V), so we ob tain A v 



5.34 x 10 - 22 Nh o . According to Shetty et al. | |201ll) and 
Ji et al| (|2012l ). we used the X factor between H2 and 
12 CO to estimate the column density of the clump. So 
Nh 2 — X x Wco[ cr n~ 2 ]i where Wco is the integrated 
12 CO intensity, W C o = J T mb dv cm' 2 Kkms' 1 . Fi- 
nally we obtained the A v is between 5 to 30 mag. 
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We selected the best SED fitting model (see Fig. 8). 
The source age is 1.81 x I0 4 yr, the source mass is 
Ah = 27.44M Q , the disk mass M disk = 4.52 x 10~ 2 M Q , 
the envelope mass M env = 1.48 x 1O 3 M , and the enve- 
lope accretion rate M env = 3.34 x 10" 3 M®/yr. N10-7 
satis fy the criteria M env /M+ > 10 _6 j/r _1 (Robitaille et 
al.J (|2007tl ). and is one stage I YSO. 



Acknowledgements This work was funded by The 
National Natural Science Foundation of China under 
grant 10778703, and partly supported by China Min- 
istry of Science and Technology under State Key Devel- 
opment Program for Basic Research (2012CB821800) 
and The National Natural Science Foundation of China 
under grant 10873025. 



cn 10" 





0.1 1 10 100 1000 

A (fim) 

Fig. 8 The SED fitting of high-mass YSO N10-7. Black 
line shows the best fit. The dashed line shows the stellar 
photosphere corresponding to the central source of the best 
fitting model, as it would look in the absence of circumstellar 
dust. The black points represent the input data at 3.6, 4.5, 
5.8, 8, 450, 850 and 1100 fim. 



4 Summary 

We studied the high-mass clump associated with the 
bubble N10 in detail based on the Spitzer data, contin- 
uum observations at 450, 850 and 1100 /im and spectral 
lines of 12 CO(3 - 2), HCO + (4 - 3), N 2 H+{A - 3) and 
CH 3 OH(7(0,7)- 6(0,6)). 

We draw some tentative results for the high-mass 
clump. Optically thick line 12 CO(3 — 2) are all blue- 
shifted with respect to the system velocity, this seems to 
indicate that the outer envelope of the high-mass clump 
is still falling toward the center. The position-velocity 
diagram of 7V2-ff + (4 — 3) indicates the cold dense core of 
the clump has not been destroyed by the star formation 
activities in it. Detection of CH 3 OH(7(0, 7) - 6(0,6)) 
suggests that a hot core has formed due to the feed- 
back of YSO N10-7. The mass of N10-7 is about 
27.44 M Q . The line wings of CH 3 OH{7{0, 7) - 6(0, 6)) 
seems to trace the outflow excited by N10-7. The ra- 
tio HCO+{4: - 3)/Y 2 iJ+(4 - 3) in the outer part of 
the clump is larger than that in the inner part of it. 
The reason may be that the CO abundance relative to 
Y 2 iJ + (4 — 3) increased in the outer part of the high- 
mass clump, more N2H + (A — 3) were converted into 



#CO+(4-3). 
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